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Abstract 

We simulate lattice QCD with dynamical u and d quarks at finite chemical potential, fxj, for 
the third component of isospin (I3), at both zero and at finite temperature. At zero temperature 
there is some \x c say, above which I3 and parity are spontaneously broken by a charged pion 
condensate. This is in qualitative agreement with the prediction of effective (chiral) Lagrangians 
which also predict [i c = m n . This transition appears to be second order, with scaling properties 
consistent with the mean-field predictions of such effective Lagrangian models. We have also 
studied the restoration of I3 symmetry at high temperature for > fi c . For fij sufficiently large, 
this finite temperature phase transition appears to be first order. As fii is decreased it becomes 
second order connecting continuously with the zero temperature transition. 
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I. INTRODUCTION 



Neutron stars are made of dense cold nuclear matter - hadronic matter at high 
baryon-number density and low temperature. Large nuclei can be considered as droplets 
of nuclear matter. The relativistic heavy-ion collisions now being observed at RHIC and 
the CERN heavy-ion program can produce hadronic matter at high temperature and finite 
baryon-number density. Nuclear matter also has a finite (negative) isospin (I3) density due 
to Coulomb interactions, and it has been suggested that at very high densities it would have 
a finite strangeness density. Hence it is of interest to study QCD at finite quark/baryon- 
number density, finite isospin density and finite strangeness density at both zero and finite 
temperature. 

Finite density is customarily studied through the introduction of a chemical potential 
for the charge of interest in the action. Introducing a finite chemical potential for quark- 
number leads to a complex fermion determinant with a real part of indefinite sign, which 
precludes use of standard simulation methods which rely on importance sampling. If we 
include a finite chemical potential, fij, for I 3 in the absence of any quark- number chemical 
potential, the fermion determinant remains non-negative, and simulations are possible. Such 
simulations can determine the QCD phase structure on one surface in the phase diagram 
for nuclear matter. One can hope that this will identify phases which will persist to finite 
baryon/quark-number density, and determine their properties. 

We have performed simulations of lattice QCD with 2 flavours of light staggered 
quarks at finite fij for both zero and finite temperatures. Preliminary results of these simu- 
lations were reported at Lattice2001 [Q. We included a small explicit /3-breaking interaction 
needed to observe spontaneous symmetry breaking on a finite lattice. In addition to allow- 
ing us to observe spontaneous breaking of I3 (and parity), this term renders the fermion 
determinant strictly positive. Our zero temperature simulations were performed on an 8 4 
lattice at an intermediate value of the coupling constant. In the limit that our symmetry- 
breaking parameter vanishes, there is some fii = /i c above which J 3 and parity are broken 
spontaneously by a charged pion condensate. This is in accord with the predictions of Son 
and Stephanov using effective (chiral) Lagrangians 0. We observe critical scaling consis- 
tent with the mean-field predictions of these effective (chiral) Lagrangians. In drawing these 



2 



conclusions, it is important to use equations of state of the forms predicted by such effective 
Lagrangians. These results are similar to those obtained in studies of the quenched version 
of this theory ||. They are also similar to what is observed for 2-colour QCD at finite 
chemical potential \i for quark- number 171] . This is not surprising since the effective 

Lagrangian analysis of 2-colour QCD at finite /i/ || is similar to that for QCD at finite \x\. 

Our finite temperature simulations were performed on 8 3 x 4 lattices. At sufficiently 
high temperature and fij > fi c , we observe the evaporation of the symmetry-breaking pion 
condensate. For fii sufficiently large, this transition is first order. As fij —>■ fi c this transition 
softens and appears to become second order. Such a transition from first to second order 
should occur at a tricritical point. Again these results are similar to what we observed for 
2-colour QCD at finite quark-number chemical potential 0. 

Section 2 gives details of the actions and their symmetries. In section 3 we present our 
zero temperature results and scaling analyses. The finite temperature results are presented 
in section 4. Discussions, conclusions and an outline of planned extensions are given in 
section 5. 

II. LATTICE ACTION AND SYMMETRIES 

The staggered fermion part of the action for lattice QCD with degenerate u and d 
quarks at a finite chemical potential \ii for isospin (7 3 ) is 

S /=E#3|ii)+mk (1) 

sites 

where is the standard staggered ^)with links in the +t direction multiplied by and 
those in the —t direction multiplied by @. When — m — 0, this action has a global 
U(2) x U(2) flavour symmetry under which 

X — ► exp[z(0 + eft). t]x 

X — ► Xexp[-z(# -e0).r] (2) 

where r = (l,r), 9 and <fi are site-independent 4-component "vectors", and e = e(x) = 
(— \Y +y+z+t . Spontaneous symmetry breaking can occur in any direction in this space. If 
we keep /i/ = and allow m/0, the symmetry is broken down to U (2)y. On the other hand 
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if we keep m = and allow \xi 7^ 0, the symmetry is broken down to U(l) x U(l) xU(l)xU (1) 
generated by 1, T3, e and er^. Finally in the general case where neither \xi nor m vanishes, 
the symmetry is reduced to U{l)y x U{l)y associated with 1 and r 3 . 

In order to predict potential symmetry breaking patterns we make several simple 
modifications of the arguments of Son and Stephanov to apply them to the staggered 
lattice action. The generic quark bilinear which creates a meson has the form 

M = xT X . (3) 

The propagator for such a meson obeys the inequality 

\(M(x)M*(0))\ < const(S(x 7 0)S t (x,0)). (4) 

Thus, meson operators M whose propagators are proportional to (S(x,0)S'(x,0)), are po- 
tential Goldstone bosons. Now we note that our Dirac operator obeys 

= [Kr 3 /i / )+m]t, (5) 

so ix er i,2X are Goldstone candidates. These are linear combinations of 7T + and tt~ creation 
operators which means that if spontaneous breaking of the remnant flavour symmetry should 
occur, one linear combination of 7r =l= will become a Goldstone boson while the orthogonal 
linear combination will develop a vacuum expectation value - a charged pion condensate. 
We note, in passing, that in the limit of massless quarks 

1-1,2^x3/17)7-1,2 = -fli^mY, (6) 

and we have 2 additional Goldstone boson candidates, x T i,2X^ an d if spontaneous symmetry 
breaking does occur we will have 2 Goldstone bosons rather than 1. 

Since, in order to observe spontaneous symmetry breaking on a finite lattice, one 
needs to add a small explicit symmetry breaking term in the direction defined by the con- 
densate, we choose to work with the fermion action 

S f = S WKn^l) + m\x + iXexT 2 x] (7) 

sites 

where the term proportional to the (small) parameter A serves this purpose. The Dirac 
operator now has the determinant 

det [pjzHi) +m + i\er 2 ] = det [A ] A + A 2 ] (8) 
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where we have defined 

A = ^i)+m. (9) 

(Note that this is a 1 x 1 matrix in the flavour space on which the rs act.) We see that adding 
this symmetry breaking term has the effect of rendering the determinant strictly positive, 
which enables us to use the hybrid molecular-dynamics (HMD) algorithm to simulate this 
theory. Note that this theory has 8 continuum flavours. We use the HMD method to take 
the required fourth root of the determinant reducing this to 2 continuum flavours. For the 
purpose of simulation, it is convenient to multiply the Dirac operator on the left by the 
matrix diag(l, — e) and on the right by the matrix diag(l,e). The transformed matrix Ai, 
has the same determinant as the original Dirac operator, and M.^ M. is block diagonal, with 
the upper and lower blocks having the same determinant. This means that we use 'noisy' 
fermions and generate Gaussian noise for both upper and lower components of Aix, but only 
keep the upper components of x after the inversion. Thus we still have only 8 flavours in 
the quadratic formulation. This is completely analogous to the odd-even lattice separation 
which prevents further species doubling in staggered lattice QCD at zero chemical potential. 

Quantities we measure include the chiral condensate, 

<W (xx), (io) 

the charged pion condensate 

i(*PlbT2^} ^i{x^2X) (11) 

and the isospin density 

Here we have included both the lattice and continuum versions of the condensates. To get 
this simple continuum form for the charged pion condensate requires absorbing a factor of 
£5 (the flavour analogue of 75) into the definition of the <i-quark field. 



III. LATTICE SIMULATIONS AT ZERO TEMPERATURE 



We have simulated Nf = 2 lattice QCD at finite fij on an 8 4 lattice at an intermediate 
coupling (3 = 6/g 2 = 5.2. This (3 was chosen since it represents an approximate lower bound 
to estimates of the finite temperature transition value for N t = 4 in the chiral limit. This 
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was used to keep finite volume effects at acceptable levels. We performed simulations at 2 
different quark masses (m = 0.025 and m = 0.05) to see that varying the mass did not affect 
the qualitative behaviour of the theory and that we understood the effects of changing the 
quark mass. 

At m — 0.025, we performed runs, each of 2000 molecular-dynamics time units in 
length, at 17 different /i/ values (0 < fii < 2) for each of A = 0.0025 and A = 0.005. Using 2 
A values, both chosen to be much less than m, enabled us to extrapolate to the A = limit, 
which is our ultimate interest. (We also ran at fii = 3.0, A = 0.005 to check saturation.) 

The charged pion condensates, i^'y&ip) from these simulations are presented in 
figure [l] as functions of fii, along with a linear extrapolation to A = 0. This extrapolation 
strongly suggests that the A — > condensate vanishes for /// < fj, c ~ 0.3 — 0.45, above 
which it is finite. This would indicate a phase transition to a phase in which /3 symmetry 
is broken spontaneously by a charged pion condensate, with an associated Goldstone mode. 

This behaviour is predicted by the effective (chiral) Lagrangian analyses of Son and 
Stephanov, which also predict that the transition should be second order with mean-field 
exponents 0. We fit our extrapolated 'data' to the critical scaling form 

^7572^) = const 0/ - ii c f m (13) 

for yU/ > /i c close to the transition. Fitting to this form over the range 0.4 < /// < 1.0, 
we find fi c = 0.394(1), f3 m = 0.230(9) and const = 1.23(2) at a 62% confidence level. This 
appears inconsistent with mean field scaling for which f3 m = ~, and closer to the tricritical 
scaling for which (3 m = |. Indeed, good fits to a tricritical scaling form over this range of \ii 
and both As can be obtained (confidence level 25%). 

However, as we have noted in our paper on the quenched theory || such fits can be 
deceptive and can be because the true scaling behaviour of these theories is best described 
by the scaling forms given by effective Lagrangian analyses. When this form of scaling is 
described in terms of hj — fi c rather than the natural scaling variables these forms imply, the 
true scaling window is very narrow. Outside this window these theories can appear to scale 
with j3 m which is half the true value when analyzed in terms of fii — fi c . We now introduce 
the scaling forms (equations of state) suggested by such effective Lagrangians, both of which 
give mean-field scaling behaviour. 
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FIG. 1: Charged pion condensate as a function of \ij for A = 0.0025, A = 0.005 and A — ► 0. The 
curves are fits of the finite A measurements to the scaling forms defined in the text. 

The form for the equation of state suggested by the lowest order tree-level analysis 
of effective Lagrangians of the non-linear sigma model type [|| is given in terms of a which 
minimizes the effective potential 
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-a /j? sin 2 (a) — bm cos(a) — b A sin(a) 



in terms of which 



^('075^2^') = b sin (a) 
(ijjip) = b cos(a) 



(14) 

(15) 
(16) 



and 

jo = 4a n sin 2 (a). (17) 
b is given in terms of \x c and a, namely 

2 

b=—aa 2 (18) 
m 

The form for the equation of state which is derived from an effective Lagrangian of 
the linear sigma model type is obtained by extracting the values of R and a which minimize 
the effective potential 

S = - R 4 — - a R 2 — -b fi 2 sin 2 (a) — cmR cos(a) — cXR sin(a) (19) 



in terms of which 



and 



c is given in terms of \i c by 



i(ip r Y5T2ip) = c R sin (a) (20) 
(tpip) = c R cos(a) (21) 

jl = 2 6 /ii? 2 sin 2 (a). (22) 
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c= b ^a + brf. (23) 



Finally, the tricritical scaling form which we use for comparison, and which does 
not yield mean-field behaviour, is expressed in terms of the value of <fi which minimizes the 
effective potential 



in terms of which, 



£ = - l - c A 3 - i (fzj - /i c ) 2 - b A (24) 



^75^> =&0+^c0 3 (25) 



and 

]o = a<f) 2 . (26) 

We fit our measurements of the charged pion condensate for < < 1 and both 
As to the linear sigma model form of equation |2(| We obtained a fit with fi c = 0.4036(5), 
a = 0.52(1), m = 0.0253(1) with x 2 /d°f — 2.2. Although worse than the tricritical fit, we 
note that it allows a reasonable fit for \ii < \i c in addition to fij > \i c which the tricritical 
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fit did not. /i c = 0.426(3) for the tricritical fit, which is close to the value expected for the 
apparent tricritical scaling produced by such sigma model scaling. One of the reasons the 
tricritical fit gives better results is that it includes a second symmetry breaking interaction, 
cubic in the order parameter, which rounds off the curve as fij approaches 1.0 allowing a 
better fit to the 'data'. Such cubic terms could be incorporated equation [T9|. However, now 
there is not one such term but several, which is why this possibility was not considered. 
Finally, as we shall see below, this fit makes good predictions for the chiral condensate and 
the isospin density. 
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FIG. 2: Chiral condensate as a function of [ij for A = 0.0025 and A = 0.005. 



In figure |2| we present the chiral condensate for the same set of simulations. The 



general characteristics of these graphs are that (ipip) remains roughly constant for fij < fi c , 
after which it drops rapidly, approaching zero at large Although this can qualitatively 
be thought of as the condensate rotating from the chiral to the isospin-breaking direction, 
it is not a simple rotation, since the magnitude of the total condensate increases up until 
saturation effects take over. This contrasts with the predictions of lowest order effective 
Lagrangians of the non-linear sigma model type where it is a simple rotation. However, in the 
case of 2-colour QCD at finite quark-number chemical potential, whose effective Lagrangian 
is structurally very similar to that for the theory at hand, the effective Lagrangian/chiral 
perturbation theory calculations have been extended to next-to-leading order [p| . Here, 
although mean field scaling survives, the rotation of the condensate is accompanied by a 
rescaling. Such behaviour is reproduced at tree level by effective Lagrangians of the linear 
sigma model variety which we use for our fits to the pion condensate. The predictions these 



fits make for the behaviour of the chiral condensate (equation |21]), have been superimposed 
on the measurements of figure [| Except for small departures at small /i/, which we attribute 
to the fact that m for the fit differs slightly from the true m = 0.025, these predictions are 
very good until the effects of saturation start to dominate. 

Figure || presents the isospin (J 3 ) density as a function of /!/. While the 2 conden- 
sates are normalized to 4 flavours for comparison with previous (fij = 0) simulations, j'q is 
normalized to 8 flavours, which is the natural normalization for staggered quarks. Qualita- 
tively, we note that j$ is close to zero for fii < fi c , rises slowly (in comparison with the pion 
condensate) up to fii ~ 1, after which it starts to rise more rapidly reaching its saturation 
value of 3 (1/2 for each of 3 colours and 2 'flavours' per site) due to fermi statistics, at 
fij ~ 2. (In fact measurements made at /// = 3.0 give a value consistent with 3). Saturation 
is clearly a finite lattice spacing effect and hence requires no further discussion. We also note 



that there is very little A dependence. The predictions of equation |22| using the parameters 
of our fit are superimposed on our 'data' and show good agreement out to jii ~ 0.8. Both 
'data' and fit are approximately linear in fij over this range. 

As also noted by Son and Stephanov, measuring j'q as a function of fii at T = and 
constant volume (since f3 is constant) yields the pressure p and energy density e as functions 
of /i/, since 

P = / fodvi (27) 
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FIG. 3: Isospin(/3) density as a function of hi for A = 0.0025 and A = 0.005. 



for hi > He and zero for \ii < /V In the region where jg ~ const (hi — He) these yield 

p = - const (hi - yU c ) 2 
1 



e — - const (/ij — h 
p _ hi - u c 
e hi + A*c 



(28) 

(29) 
(30) 
(31) 



The last of these equations is a form of the equation of state for this system in the neigh- 
bourhood of He- Clearly we could extend each of these expressions beyond the scaling region 
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by interpolating the 'data' of figure ^| and performing the relevant integrals analytically or 
numerically. 

SU(3) N f =2 £ = 5.2 m=0.05 8 4 lattice 
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FIG. 4: Pion condensates as functions of fxj for A = 0.005, A = 0.01 and a linear extrapolation to 
A = 0.0. The lines are the tricritical fit described in the text. 

We performed similar simulations at the same coupling j3 = 5.2 and mass m = 0.05 at 
A = 0.005 and A = 0.01. Here we concentrated on the neighbourhood of the phase transition 
and used more closely spaced (in //j) points with somewhat lower statistics. For fij > fi c we 
find acceptable fits of the pion condensate to to both the non-linear sigma model scaling form 
and to the tricritical scaling form (34% and 40% respectively). The fit to nondinear sigma 
model effective Lagrangian scaling (equation [15D enables one to extend this to low A fit 
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of the 'data' over the complete range over which we have measurements at both A values - 
0.4 < m < 0.7 — yields /m c = 0.569(2), a = 0.0868(9) and m = 0.0534(3), with x 2 /dof = 1.5 
for the non-linear sigma model form. Even though this fit is worse than those restricted to 
Hi ~ He-, we feel that the ability to fit fij < fi c in addition to fij > /i c makes the argument 
for this form of fit more compelling. Not only do we see qualitative consistency with the 
smaller mass results, but our measured values of /i c are consistent with the expectation that 
/i c (m = 0.05) = \/2[i c {m = 0.025) which would be true if, indeed, fi c = m^, from PCAC. 
This 'data' for the pion condensate with the scaling fits superimposed is plotted in figure |j. 
We note that the fits appear to have validity beyond the range of \ii over which the fits were 
performed. 



IV. SIMULATIONS AT FINITE TEMPERATURE 



We have performed simulations of QCD at finite /i/ and finite temperature on an 
8 3 x 4 lattice, with m = 0.05. Most of these simulations were performed at A = 0.005 and 
A = 0.01, i.e. A << m, with the objective of obtaining information about the A = limit. 
The goal of these simulations is to map out the region of the (f3, /i/) and hence the (T, /xj) 
plane where I 3 is spontaneously broken by a charged pion condensate, and determine the 
nature of the phase transitions which demarcate its boundaries. 

The first of these simulations was performed at a fixed, large (but well below satu- 
ration) value of /!/. The value chosen was /ij = 0.8. At (3 low enough to approximate zero 
temperature on an N t = 4 lattice, the system is in the phase where I3 is spontaneously 
broken by a (large) pion condensate. As (3 is increased we eventually reach a value (3 = (3 C at 
which this condensate evaporates. For (3 > (3 C we are in the phase where the pion condensate 
vanishes for A — ► and I3 is unbroken. A single A value, A = 0.005 was used for these runs. 

Figure |5| shows the f3 dependence of the pion condensate for these simulations. We 
see that for (3 < 5.2 the condensate is large. Between {3 = 5.2 and (3 = 5.3 the condensate 
drops by an order of magnitude, and is small enough for (3 > 5.3 that we are safe to assume 
that it would vanish in the A — > limit. This drop is so precipitous that we suspect that 
it is first order, although we really need a larger lattice to confirm this. 

We have also measured the Thermal Wilson Line (Polyakov Loop) during these runs. 
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FIG. 5: Pion condensate as a function of (3 for m = 0.05, A = 0.005, and = 0.8 on an 8 3 x 4 
lattice. 

These measurements are shown in figure |6|. For f3 < 5.2, the Wilson line is small indi- 
cating confinement. For (3 > 5.3, the Wilson line becomes large indicating deconfinement. 
The jump between (3 = 5.2 and (3 = 5.3 is again great enough to suggest a first order 
transition. This behaviour of the Wilson Line indicates that this is the temperature-driven 
deconfinement transition. 

We have also run simulations on an 8 3 x 4 lattice with f3 = 4.0 which gives us the low 
temperature behaviour. We chose m = 0.05 again and ran at A = 0.005 and A = 0.01 for 
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FIG. 6: Wilson linear a function of f3 for m = 0.05, A = 0.005, and [ii = 0.8 on an 8 3 x 4 lattice. 

0.20 < hi < 2.0 which covers both the transition from the I3 symmetric phase to the phase 
where J3 is spontaneously broken, and the approach to saturation. Our measurements of 
the pion condensate are shown in figure [7[ Not surprisingly this graph resembles those we 
obtained on an 8 4 lattice since at (3 = 4.0, this lattice is essentially at zero temperature. 

We fit the scaling behaviour of these measurements to the non-linear sigma model 



scaling form of equation [15] for both A values and 0.2 < fij < 0.8. The fit yielded \i c = 
0.519(1), a = 0.1400(4) with the quark mass fixed at m = 0.05 at a confidence level of 41%, 
which is very good. These fits are shown in figure |?]. Again an acceptable tricritical fit was 
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FIG. 7: Pion condensates as functions of fxj for A = 0.005, A = 0.01 and a linear extrapolation to 
A = 0.0 on an 8 3 x 4 lattice at (3 = 4.0. The lines are the tricritical fit described in the text. 

possible (confidence level 48%), but only for 0.55 < /ij < 0.8, yielding the expected larger 
estimate for /i c . 

In figure |8] we present the isospin density from these runs, j'q rises from zero at 
Hi ~ /i c . Again there is little A dependence. We have superimposed the form predicted from 
the fit to pion condensate using equation [T7| on these plots. These curves are in reasonable 
agreement with the measured values up to /i/ ~ 0.7. This indicates that the scaling window 
for j'q is slightly less than that for the pion condensate. However, it is clear that the linear 
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FIG. 8: Isospin(/3) density as a function of fij for m = 0.05, A = 0.005 and A = 0.01 on an 8 3 x 4 
lattice at j3 = 4.0. 

or near-linear increase of this quantity with \xj continues beyond the point where the 'data' 
and curves diverge. This is born out by the fitting j\ to the form 



jo = const {fx i - ll c ) Pi . 



(32) 



Fitting the A = 0.005 'data' over the range 0.5 < fij < 1.1 gives li c = 0.467(13), & = 0.94(6) 
and const = 1.22(3) at a confidence level of 18% while fitting the A = 0.01 'data' yields 
ii c = 0.382(9), Pi = 1.09(3), const = 1.20(1) at a confidence level of 85%. These results 
are in excellent agreement with the effective Lagrangian prediction f3j = 1. This graph also 
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indicates a crossover to a more rapid increase at \xj ~ 1.5. 
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FIG. 9: Pion condensates as functions of fij for A = 0.005, A = 0.01 and a linear extrapolation to 
A = 0.0 on an 8 3 x 4 lattice, at (3 = 5.0. The lines are the tricritical fit described in the text. 

Finally we have performed m = 0.05 simulations on an 8 3 x 4 lattice at (3 = 5.0, which 
lies below f3 c at high fij, while being large enough that the effects of finite temperature should 
be apparent. Here again we ran with A = 0.005, 0.01. Here we used higher statistics (2000 
time-units per 'run') but at fewer fij values. Our measurements of the charged pion con- 
densate are presented in figure ^. The transition from low to high values of this condensate 
shows no sign of a first order transition. The linear extrapolation to A = gives values 
close to zero at low fij, rising rapidly from zero above some \i c ~ 0.5. Here we find good 
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fits of the 'data' for both A values and 0.2 < fij < 1.0 to the non- linear sigma model form 
with fi c = 0.5521(5), a = 0.1040(4), m = 0.0523(1) at a confidence level of 31% and to the 
linear sigma model form with /i c = 0.5513(5), a = 2.46(1), b = 0.0823(4), m = 0.0518(1) at 
a confidence level of 48%. This second fit in included in figure |^. A tricritical fit is possible 
for 0.55 < /ij < 0.8 and has a confidence level of 28%. 

Thus the line of phase transitions which bound the region in the (/!/, T) plane within 
which isospin (^3) is broken spontaneously by a charged pion condensate, is second order for 
low temperatures and becomes first order at high [ij. The second order segment of this line 
appears to have mean-field critical exponents. 



V. DISCUSSION AND CONCLUSIONS 



We have simulated QCD with 2 quark flavours (u,d) at a finite chemical potential (fij) 
for isospin, I 3 . At zero temperature and intermediate coupling (/3 = 6/g 2 = 5.2) we found 
strong evidence for a second order transition to a phase in which I 3 is spontaneously broken 
by a charged pion condensate which also breaks parity, at /!/ = /i c . The observed behaviour 
is what is predicted by effective Lagrangian methods for fif appreciably less than the value at 
which saturation, a lattice artifact, takes over. These effective Lagrangian analyses predict 
that /i c = m^. Since we have not measured on these small lattices, all we have checked 
is that /i c oc \fm for the 2 quark masses which we use (m = 0.025,m = 0.05), and found that 
this is true within the uncertainties of our measurements. The critical scaling appears to 
be well described by an equation of state suggested by these effective Lagrangian analyses, 
which means that the critical point has mean-field critical exponents. However, tricritical 
behaviour cannot be completely excluded. Such behaviour was discussed in detail in our 
paper on the quenched theory [[|. 

We also measured the isospin density (j'q) and found that it increases from zero, at or 
near /i c . The scaling behaviour appears to be linear close to predicted by effective 

Lagrangian analyses M and well described by the predictions given by the fits to the pion 
condensate, within the scaling window. For larger fij values it starts to increase considerably 
faster than linear. This is in qualitative agreement with the expectation that J3 density 
should increase as /zf at large [ij. We have not been able to determine if our observations 
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are consistent with this /xf increase because of the effects of saturation which cause the 
isospin density to approach 3 at high fij. Here we have indicated how the measurement of 
jo enables one to obtain the pressure (p) and the energy density (e) as functions of ///, and 
given explicit expressions for these quantities in the scaling regime. 

The chiral condensate remains approximately constant for \x\ < fi c . Above \x c it starts 
to fall approaching zero for large /i. Again this is in agreement with expectations, and the 
predictions from the fits to the pion condensate. However, the expectation from lowest-order 
effective Lagrangian tree-level analyses, that the chiral condensate simply rotates into the 
direction of the charged pion condensate, is not true. However, in closely related 2-colour 
QCD at finite quark- number chemical potential, chiral perturbation theory calculations 
through next-to-leading order show that, while scaling remains mean field, the condensate 



does not merely rotate, but also rescales |TT[. Since the structure of chiral perturbation 
theory (effective Lagrangians) for the two theories is so similar, we expect a similar result 
for QCD at finite /ij. 

We have performed simulations at finite temperature (T) in addition to finite /ij. In 
particular we have heated the system at fixed \xi = 0.8 > n c {T = 0) (in lattice units) by 
increasing (3. On our N t = 4 lattice there is some f3 = (5 C (5.2 < (3 C < 5.3) at which the 
charged pion condensate evaporates. This transition appears to be first order. Such first 
order behaviour was predicted for /i/ large enough by Son and Stephanov J|] who argued 
that at high fij the fermions would effectively decouple, and the phase transition would be 
that for pure glue, which is known to be first order. We note from our observations that 
the situation is not quite this simple. The pure glue transition on an N t = 4 lattice occurs 



at (3 ~ 5.7 JL2], |i3[ . Since our observations place (3 C somewhat lower than this, and close 
to the value of the /// = finite temperature transition, the quarks are having an effect. 
We also note that Son and Stephanov suggest that the first order deconfinement transition 
at high temperature is distinct from the -Z3 symmetry restoring transition. Our evidence 
is that these 2 transitions are coincident. This means that the second order segment of 
the phase boundary can be considered as driven purely by the chemical potential, which 
makes their argument for 0(2) universality less compelling. We note that, in our paper on 
2-colour QCD at finite fi and T 0, we present an alternative argument for the first-order 
finite-temperature transition which that theory exhibits for large /1. 
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In addition, we have simulated our N t = 4 system at fixed j3, varying /xj. In particular 
we have performed simulations at /3 — 4.0 which is at near-zero temperature, and (3 = 5.0 
where the system is clearly at a finite temperature. Both of these simulations showed a 
second order transition. Again the scaling was well described by the scaling forms suggested 
by effective Lagrangian analyses which indicates that they have mean-field critical exponents. 

We have noted throughout this paper the similarity between 3-colour QCD at finite /xj 
and 2-colour QCD at finite quark-number chemical potential, /x. The correspondence is seen 
by identifying |/xj with /x, the charged pion condensate with the diquark condensate and the 
isospin density with the quark number density. This similarity is seen both in simulations and 
in effective Lagrangian analyses. The expected position of the zero temperature transition 
is the same /xj = m n (/x = \m-^). Its nature — second order with mean-field exponents 
- is the same as predicted |L1| and observed || [7| for 2-colour QCD. In both systems the 
spontaneous symmetry breaking is in the Goldstone mode (superfluid). At finite temperature 
the condensate evaporates at a transition which is first order for /// (/x) large enough. This 
line of first order transitions softens to second order and the line of second order transitions 
connects to the zero temperature transition. Thus, we can use 2-colour QCD results at finite 
/x as a guide as to QCD at finite /xj. 

We are extending these simulations to a larger lattice (12 3 x 24) and weaker coupling 
where we hope to observe the expected mean-field transition more clearly distinguished from 
the tricritical alternative and rule out the 0(2) alternative. This lattice will also enable us 
to measure the spectrum of Goldstone and pseudo- Goldstone excitations as functions of /X/. 
More extensive spectrum analyses at fij = will give us a more definitive scale for these 
phenomena, in addition to a value for with which to compare /x c . Configurations will 
be stored so that we can make other spectroscopy measurements at finite /xj. We will also 
extend the finite temperature simulations to a 12 3 x 6 lattice since it is difficult to determine 
the order of the continuum transitions from 8 3 x 4 lattices. In addition we will study the 
instantons at large /xj since it is believed that instantons and their interactions have a 
relatively simple structure at large isospin density, analogous to what has been predicted for 
large quark- number density \T%\. 

We are also starting lattice QCD simulations including both a chemical potential \xi 
for isospin and a chemical potential \x s for strangeness. Here effective Lagrangian analyses 
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have indicated that there is a competition between the formation of pion and kaon conden- 
sates as \ii and \i s are varied independently and that the boundary between the region with 



a pion condensate and that with a kaon condensate is a line of first order transitions [ R ] . 

Although our simulations have been limited to zero baryon number density, it is 
interesting to know how much of this analysis is relevant to systems with finite baryon 
number density and isospin density. If it does have relevance, charged pion condensates could 
contribute to the equation-of-state of nuclear matter and thus be important in understanding 
the physics of neutron stars and perhaps large nuclei. 
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